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ABSTRACT 

The young and compact planetary nebula IC 4997 is known to exhibit very broad 
wings with a width exceeding 5000 km s~^ around Ha. We propose that the broad wings 
are formed through Rayleigh-Raman scattering involving atomic hydrogen, by which 
Ly/3 photons with a velocity width of a few 10^ km s~^ are converted to optical photons 
and fill the Ha broad wing region. The conversion efficiency reaches 0.6 near the line 
center where the scattering optical depth is much larger than 1 and rapidly decreases 
in the far wings. Assuming that close to the central star there exists an unresolved 
inner compact core of high density, uh ~ 

109-10 

cm ^, we use the photoionization code 
'CLOUDY' to show that sufficient Ly/3 photons for scattering are produced. Using a 
top-hat incident profile for the Ly/3 flux and a scattering region with a H I column 
density Nhi = 2 x 10^*^ cm~^ and a substantial covering factor, we perform a profile 
fitting analysis to obtain a satisfactory fit to the observed flux. We briefly discuss the 
astrophysical implications of the Rayleigh-Raman processes in planetary nebulae and 
other emission objects. 

Subject headings: planetary nebulae — planetary nebulae : individual IC 4997 — ra- 
diative transfer — scattering — profile 

1. Introduction 

With the advent of the Hubble Space Telescope and large optical telescopes, planetary nebulae 
(PNe) came to be known to possess complex structures including multiple shells and collimated 
outflows, of which the origins still remain puzzling (e.g. Gurzadyan 1997). In addition, PNe 
display various nebular morphology ranging from elliptical to spherical, bipolar (or of butterfly 
shape), point-symmetrical and irregular. However, the physical explanation for the morphological 
diversity is not found (e.g. Corradi & Schwarz 1995). 
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Noting that symbiotic stars, widely known to be binary systems of a cool giant and a hot white 
dwarf or a main sequence star, exhibit bipolar emission nebulae, Soker (1998) advocated the view 
that the bipolar morphology in planetary nebulae is attributed to the binarity in the central star 
system. Evolutionary links between symbiotic stars and PNe have been also proposed by a number 
of researchers (e.g. Corradi 1995). Currently one of the most controversial issues is the existence of 
an accretion disk in symbiotic stars and bipolar PNe. According to the recent SPH computations 
by Mastrodemos &; Morris (1998), in a detached binary system consisting of a mass losing giant 
and a companion star, the slow and dusty stellar wind from the giant may induce the formation of 
an accretion disk around the companion star (see also Theuns et al. 1996). 

Circumstantial evidence supporting this scenario was provided by Lee &; Park (1999), who 
proposed that a disk type emission model can naturally explain the main characteristics of the 
profiles and the polarized fiuxes of the Raman-scattered 6830 and 7088 A features in the symbiotic 
star RR Tel. The Raman scattering process by atomic hydrogen has been introduced by Schmid 
(1989), who identified the 6830 and 7088 A bands in symbiotic stars as the Raman scattered O VI 
1032, 1038 doublet. A thick H I component around the giant companion in a symbiotic star is 
believed to act as the Raman-scatterer which converts the UV continuum around Ly/3 into the 
optical continuum. 

The young and compact PN IC 4997 is particularly interesting in that it also displays a hint 
of bipolar morphology, as is revealed in the radio image presented by Miranda & Torrcllcs (1998). 
Miranda et al. (1996) presented spectroscopy of IC 4997 which shows very broad wings around Ha 
with FWHM 5375 km s"^. IC 4997 is not the only PN showing the broad Ha wings. The bipolar 
planetary nebula M 2-9 is also known to exhibit extremely broad wings with width ~ 10^ km 
in Ha (Balick 1989, see also Torres-Peimbert Sz Arrieta 1998). 

Patriarchi and Perinotto (1991) used the lUE data of 159 planetary nebulae to find that about 
60 per cent of central stars have a stellar wind with a measured edge velocity ranging 600 to 
3500 km s~^. Feibelman (1982) presented lUE spectroscopy of IC 4997 and found no evidence of 
P-Cygni profile in C IV 1550 doublet, which strongly implies the absence of strong stellar winds 
around the central star. This makes it difficult to attribute the broad wing around Ha in IC 4997 
to the kinematics associated with the fast stellar wind. 

Electron scattering has been often invoked to explain the wing formation in luminous blue 
variables (Bernat & Lambert 1978). However, it does not appear to be the case in IC 4997 since 
the Thomson scattering cross section is independent of wavelength and therefore broad wings 
should be present around other emission lines. Furthermore, in a planetary nebula, it is difficult 
to incorporate an electron scattering component with the electron temperature Tg > 10^ K and a 
sufficient scattering optical depth. 

An interesting observation of IC 4997 is provided by Altschuler et al. (1986), who reported 
the detection of H I in IC 4997 by observing the 21 cm absorption dip (see also Schneider et al. 
1987). With a large amount of atomic hydrogen, the broad wings around Ha can be formed by 
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Rayleigh-Raman scattering, where the incident radiation in the vicinity of Ly/3 converted to the 

radiation around Ha with a width broadened by a factor of \Ha/^LyP = 6.4 due to the incoherency 
of the scattering process. In this Letter, we investigate the broad Ha wing formation through 
the Rayleigh-Raman scattering process and discuss the imphcations on the astrophysics of bipolar 
planetary nebulae. 

2. Model 

2.1. Wing Formation from the Rayleigh-Raman Process 

In this work, we adopt a simple model to investigate the wing formation through the Rayleigh- 
Raman process. We assume that the scattering region is in the form of a finite slab with the column 
density Nhi with the solid angle il^ with respect to the radiation source, which coincides with the 
central star of the planetary nebula. We discuss the Ly/? line radiation source in the next subsection 
and now we assume that the scattering region is illuminated by a far UV radiation source with 
intensity I\. 

A UV photon around Ly/? is incident on a hydrogen atom in the ground state Is, which 
subsequently de-excites either to the 2s state re-emitting an optical photon around Ha (Raman 
scattering) or to the Is state resulting in an outgoing UV photon with the same frequency (Rayleigh 
scattering). Depending on the scattering optical depth, a UV photon may suffer a number of 
Rayleigh scatterings followed by a Raman scattering and escape the scattering region, because the 
2s state is usually populated negligibly compared with the ground state. 

The incident wavelength and the wavelength Aq of the Raman-scattered radiation are related 
by the energy conservation, A~^ = A" ^-t-A^^^,, which yields the conversion of the wavelength interval 

AA„/A„ = (A„/Ai)(AAi/Ai). (1) 

Because of the factor \o/\, UV photons with a width AV^ centered at Ly/3 are spread around Ha 
with a broadened width AV^ = 6.4AVi. 

The number of UV photons incident on the scattering region per unit time per unit wavelength 
interval is given by Ix.{T^)TrR^^s/{hc/Xi). We introduce the parameter CR{Xi), which is the con- 
version factor from UV to optical through the Rayleigh-Raman scattering process. The conversion 
factor Cii{Xi) depends sensitively on the wavelength and the column density as well as the detailed 
geometry of the scattering region, and can be computed using a Monte Carlo technique (see Lee Sz 
Yun 1998). 

The number of Raman-scattered photons coming out from the region per unit time per unit 
wavelength interval is given by 



NHa = [hi{T.)nRl^s/{hc/\)]CR{\)\y\l 



(2) 
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Therefore, the Raman-scattered flux is 

F««-(A) = NHa{hc/\)/T-' 

= F^^{X,)xC,ff (3) 

where F^^{Xi) = [Ix^7^Rl/r^] is the flux that would be observed from the radiation source with 
a size i?^, located at a distance r and Cgff = [^gCRiXi) \f / X^] is the effective Raman conversion 
efficiency. 



2.2. Conversion Efficiency 

In Fig. 1, we present the Raman conversion efficiency C/j(A) for H I column densities N^j = 
10^^,10^°, and 10^^ cm^^^ using the Monte Carlo code developed by Lee & Yun (1998). In the 
optically thin regime, the single scattering approximation holds so that 

Cr{X) = Ttot{(^Ram/(^tot)- (4) 

In the optically thick limit, the conversion rate is investigated by Lee & Lee (1997), who proposed 
an empirical relation 

Cr{X) = ^r./(n)/^[(l - r,)l3{n) + T,]f{n), (5) 

n n 

where = CFRam/f^tot- Here, (5{n) ~ \e~^ is the escape probability from the nth scattering site, 
and /(n) is the photon number flux scattered no less than n times which is obtained recursively by 

/(n + l) = (l-r,)[l-/3(n)]/(n). (6) 

A direct substitution with rv = 0.11 near the line center yields Cr{X) = 0.6, which is in excellent 
agreement with the result shown in Fig. 1. 



2.3. Photoionization Structures of the Emission Nebula 

In this subsection we discuss the Ly/5 flux that may enter the scattering region. Hyung et 
al. (1994) performed a photoionization computation to conclude that the observed emission lines 
are reproduced approximately using a model consisting of a thin shell of density n ~ 10^ cm *^ 
surrounded by a much larger shell with a much lower density ~ 10^ cm""^. They also noted that the 
N III] A1754/A1749 ratio may indicate the existence of a much denser region with n > 10^ cm~^ 
(e.g. see Czyzak et al. 1986). 

We use the photoionization code 'CLOUDY 90.05' developed by Ferland (1996) to obtain the 
Ly/3 flux that is expected around a central star of a PN. We fix the temperature and radius of the 
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central star to be = 6 x 10^ K, i?* = Rq and compute the hydrogen emission line fluxes for 

densities n = 10^ — lO^'^ cm~^. According to Hyung et al. (1994) the observed Ha luminosity is 4 x 
-j^q35 gj.g g-i assuming the distance d = 2.4 kpc to IC 4997 and the interstellar extinction parameter 
C = log I(H/3)/F(H/3) = 0.8 where the observed flux F(H/9) in H/3 = 2.95 xlQ-^^ ergs s"^ cm'^. 

In Table 1, we show the result. This shows that the observed Ha flux is obtained from a 
region of density n ~ 10^ cm for which we get a much smaller Ly/? flux. However, when 
n ~ 10^~^° cm~^, a much larger Ly/3 flux is obtained. Therefore, if there exists an emission region 
of density n ~ 10^~^° cm~^, we may have a sufflcient number of Ly/3 photons to fill the Ha wings. 
Furthermore, this region is plausibly located much closer to the central star than other emission 
regions are. Hence assuming that the distance is ~ 0.1 AU, we may expect that there is a high 
velocity dispersion of order 100 km in the region. 

3. Result 

In Fig. 2, by the dotted lines we show the Ha wings expected from the Rayleigh-Raman 
processes. For the incident Ly/3 flux we use L^yp = 6 x 10^^ erg s~^, which is expected to be 
obtained from an unresolved core of high density n ~ 10^~^° cm~^ near the central star. We tested 
two types of proflles for the incident Ly/3 flux, that is, a top-hat proflle with a width of 300 km s"-*^ 
and a Gaussian with the same velocity width. We present the results for the best-fltting column 
density of the scattering region which is Nhi = 2 x 10^*^ cm~^ for the incident top-hat profile and 
Nhi = 4 X 10^*^ cm^^ for the Gaussian profile. In this work, we concentrate only on the wing parts 
farther than 300 km s^^ (AA > 6.5 A around Ha) from the line center. 

Because the top-hat profile for the incident Ly/3 provides an excellent fit to the observed Ha 
wing profile, we expect that other types of input profiles can also produce a good fit except for some 
extreme cases. This is illustrated by another good fit using a Gaussian profile shown in Fig. 2(b). 
However, we realize the much higher column density, i.e. by a factor of 2, is required in the latter 
case. 

In this work, the incident Ly/3 flux was computed with a highly simplifled assumption. An 
introduction of a velocity fleld in the emission region may alter drastically both the strength and 
the profile of the emergent line flux. Furthermore the scattering region is complex and consists of 
a number of components differing in column density and geometrical shape. Therefore, a better fit 
can be assured from composite scattering region models and more complicated radiative transfer, 
which are deferred to a future investigation. Since the extreme wings are formed by the scatterers 
located in the equatorial region, we predict stronger polarization in these regions than in the near 
center regions that are formed by scatterers distributed nearly spherically. 
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4. Discussion 

The exact location of the scattering region with Nhi = 10^*^ cm^^ is not well-constrained 
from this analysis. If the bipolar morphology has originated from a binary central star system as 
in symbiotic stars advocated by Soker (1998), the binary companion of IC 4997 may provide a 

scattering site in the form of an extended stellar atmosphere or a slow and probably dusty stellar 
wind. This situation is often met in symbiotic stars, where a slow and thick wind from a giant 
companion may act as a Raman scattering site. In this picture, the column density of ~ 10^*^ cm~^ 
is more plausibly attributed to the slow stellar wind. 

In the bipolar planetary nebula M2-9, Balick (1989) observed the spectrum similar to that 
of the symbiotic star RR Tel. M2-9 might be the case of a binary central star system. Since 
there exists a thick circumstellar region in M2-9, the extremely wide Ha wings with a width 
11000 km are likely to be formed through the Raman process in a scattering region of column 
density Nhi ~ 10^^ cm~^. 

Broad Ha wings are found in a number of post-AGB stars according to a recent report by 
Van de Steene et al. (1999), who excluded several theoretical possibilities as the wing formation 
mechanism including Thomson scattering and Rayleigh scattering. Considering that in post-AGB 
stars much H I can be found in the stellar wind region and/or in a dense torus surrounding the 
central star, it is very convincing that the broad Ha wings in these systems are also attributable to 
Rayleigh-Raman scattering. So far the Raman-scattering is unique to symbiotic stars and the only 
exception appears to be the young planetary nebula NGC 7027, from which a Raman-scattered 
He II line has been identified (Pequignot et al. 1997). Prom this study broad Ha wings are expected 
to be found in a hot and dense H II region surrounded by an H I region, where the conditions may 
be met in objects including symbiotic stars, supernova remnants, and galactic supershells. 

Since the Rayleigh-Raman scattering provides a strong constraint on the amount of H I in the 
nebula, it can be an important tool in studying the mass loss process in the late stage of AGB 
evolution. Through high quality spectra around Balmer emission lines, the physical parameters of 
the central star can be estimated, which is complementary to the information obtained by other 
methods such as the Zanstra method. We also note that spectropolarimetry can provide more 
detailed structural information as mentioned in section 3. 
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Fig. 1. — Raman conversion efficiency Cr{X) for H I column densities Nhi = 10^^ (solid line), 
10^° (dotted line), and 10^^ cm~^ (long dashed line) obtained using a Monte Carlo technique. The 
horizontal (lower) axis is the outgoing wavelength around Ha and in the upper part is shown the 
incident wavelength around Ly/?. 
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Fig. 2. — Observed flux around Ha and the wings produced numerically in this work. The flux, 
represented by the solid line, was observed on August 31, 1991 UT, and is corrected for atmospheric 
and interstellar extinctions [C = 0.8). The vertical axis represents the specific luminosity with the 
assumed distance 2.4 kpc to IC 4997. By the dotted lines we show the wings produced via Rayleigh- 
Raman scattering. Wc assumed that the total incident Ly/3 flux L^yQ = 6 x lO''^"'' erg s^^. In panel 
(a), the incident Ly/3 has a top-hat profile with a width 300 km s^^ and the H I column density of 
the scattering region Nhi = 2 x 10^^ cm~^. In panel (b), the incident profile is a Gaussian with 
the same width and the scattering region has Nhi = 4 x 10^° cm"^. We also illustrate the flux 
multiplied by 1/40 to show the whole Ha profile. 
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Table 1. HI emission line fluxes from a plasma photoionized by a PN central star with 







= 6 X 10^ 


K, and i?* 


= Rq- 




logn^ 






log FLy^ 


log FHa 


log Fh/3 


10 


36.849 


35.687 


35.154 


35.382 


34.988 


9 


36.944 


35.233 


34.544 


35.874 


35.451 


8 


36.957 


34.025 


33.135 


35.957 


34.991 


7 


36.961 


33.268 


33.043 


35.890 


35.318 


6 


36.953 


33.146 


33.049 


35.875 


35.380 



All quantities are measured in cgs units. 



Note. — The photoionization code 'CLOUDY 90.05' (Fer- 
land 1996) was used to obtain the result. 



